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Abstract 

Practical methods for quantitative analysis of radial and angular coordinates 
of leafy organs of vascular plants are presented and applied to published 
phyllotactic patterns of various real systems from young leaves on a shoot 
tip to florets on a flower head. The constancy of divergence angle is borne out 
with accuracy of less than a degree. It is shown that apparent fluctuations in 
divergence angle are in large part systematic variations caused by the invalid 
assumption of a fixed center and/or by secondary deformations, while random 
fiuctuations are of minor importance. 

Keywords: phyllotaxis; logarithmic spiral; parastichy lattice; Helianthus; 
Asteraceae 

1. Introduction 

It has long since been recognized that divergence angles between succes- 
sive leafy organs of vascular plants are accurately regulated at one of special 
constant values. Deviations from the the constant angle are normally so small 
that this botanical phenomenon, phyllotaxis, should rather be regarded as 
a genuine subject of exact science. Besides the angular regularity, radial 
coordinates also exhibit mathematical regularity. 

In a polar coordinate system, position of the n-th leaf is specified by the 
radial and angular coordinates (r„,6'„). The angular regularity is expressed 
by the equation 

On = nd, (1) 

where d is a divergence angle. A spiral pattern is made when the radial 
component r„ is a monotonic function of n, i.e., preserves the order of the 
leaf index n. Particularly important is a logarithmic spiral given by 

Tn = a", (2) 
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where a is a constant, called the plastochron ratio. Leaf primordia at a shoot 
tip appear to be arranged on a logarithmic spiral. Eq. is a solution of the 
differential equation of exponential growth, 



dn 



log r„ = log a. 



(3) 



Thus, the logarithmic spiral is characterized by the constant growth rate 
log a. On the other side, a power-law spiral given by 

(4) 



n 



has b een oft e n use d as a m a thema t ical model of packe d seeds on a sunflower 
head flVoeell fll979h : iRidlevI fll982[ ): iRivier et all f ll984h l The equation (g]) is 
expected when all seeds have equal areas, i.e., for the area per seed. 



_(7rr„) = const. 



(5) 



In mathematical studies, the relations (II]), ([2]) and are often accepted 
without inquiring their empirical basis. From a physical standpoint, the 
validity of the mathematical relations has to be checked against experiments 
at all events. While embarking on quantitative assessment, however, we 
are confronted with a serious methodological problem. First of all, a center 
of the pattern, the origin of the polar coordinate system, must be located. 
Unfortunately, it is not always possible to locate a fixed center properly and 
objectively. The uncertainty of the central position becomes a fundamental 
obstacle to assessing the empirical relations in a quantitative manner. 

Quantitative analysis of phyllo taxis has not been made until recently. 
Maksymowych and EricksonI (119771 ) used a trigonometric method for evalu- 
ating divergence angle from positions of three successive leaves. Their three- 
point meth od is based on t he ex pone ntial gro wth ((2]), and it has been de- 
veloped by Meicenheimer ( 1986 ) and HottonI (j2003i) . On the other side, 
Rutishausen ( 1l998l ) has estimated the plastochron ratio a from the geo- 
metrical mean evaluated for selected leaves. The exponen- 
tial growth i s veri fied by the observation that the mean is nearly constant. 
Ryan et al.l ( ll99l[ ) has assessed the validity of (jl]) for real sunflowers by in- 
direct means of evaluating the Euclidean distance between successive seeds, 

= ''^n + ^n+i ~ 2r„r„+i cos(6'„ — 9n+i). This method has a merit of being 
independent of the choice of the center position. T hey have attempted pa- 
rameter fittings with various functional forms of r„. iMatkowski et al.l (Il998[ ) 



2 



has investigated computatio nal methods b y regarding the center of gravity 
as the center of the pattern. iHottonI (120031 ) has proposed to use the minimal 
variation center that minimizes the standard deviation of divergence angles. 

The main aim of this paper is to present practical methods for evalu- 
ating radial and angular coordinates of phyllotactic units without assum- 
ing any specific functional form of r„ nor the existence of a fixed center 
of the pattern. In Sec. 12. 1^ a floating center of divergence is defined by 
positions of four consecutive leaves. In Sec. 12. 2[ a floating center of paras- 
tichy is defined by four nearby points forming a unit cell of a parastichy 
lattice. For this purpose, a systematic index system for leaves of a mul- 
tijugate pattern is proposed in Sec. 12.2.21 In Sec. [31 effects of uniform 
deformation are considered. In Sec. IH presented methods are applied to 
real representative systems, which a r e arbitrarily chosen from previou s work s 
(iMaksymowvch and EricksonI ( 119771 ) : lRyan et al.l (Il99l[ ) : iRutishauserl ( 119981 ): 
Hotton et al.l (l2006l )). In [Append ix Al mathematical characteristics of the 
logarithmic spiral pattern are examined carefully. To assist in interpreting 
the results in Sec. [U general mathematical formulas for the plastochron ra- 



tio a of systems with orthogonal parastichies are derived. In [Appendix B 



practical method for numbering packed florets of a high phyllotaxis pattern is 
presented. The method is based on a numjjer th eoretic algorithm. To analyze 
phyllotaxis of Asteraceae, [Hotton et al.[ ( 12006[ ) have introduced a primordia 
front based on the assumption that the pattern has a definite center and the 
radius r„ from the center is a monotonic function of n. In Appendix B a 



similar concept, a front circle on a parastichy lattice, is introduced without 
referring to the center. 



2. Floating center 

2.1. Floating center of divergence 

For the purpose of this paper, it is suffice to investigate ph yllotactic pat- 
terns projected on a plane. For three-dimensional analysis, see [Hotton et al. 
( 12006[ ) and references therein. An effect due to a slight inclination of the 
normal axis is discussed in Sec. [31 

Consider four consecutive points Pi, P2, P3 and P4 in the same geometric 
plane. Given the first three points Pi, P2 and P3, the point O2 satisfying 



IP1O2P2 = IP2O2P3 



(6) 
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Figure 1: Given four consecutive points, Pi, P2, P3 and P4, a center of divergence O23 
is defined such that /.P1O23P2 = /-P2O23P3 = /-P3 023Pa- A dashed curve C2 is a locus 
of the point O2 satisfying ZP1O2-P2 = ^^202^3- Similarfy, a dashed curve C3 is drawn 
by the point O3 satisfying ZP2O3P3 — IP3O3P4. The point where C2 and C3 intersect is 
O23, the floating center of divergence. 



defines a curve C2 (Fig. [T]). When P2P1 = -P2-P3, or the triangle P1P2P3 is 
isosceles, C2 bisects the vertex angle IP1P2P3. Similarly, for the three points 
P2, P3 and P4, we obtain a curve C3 traced by the point O3 satisfying 

ZP2O3P3 = ZP3O3P4. (7) 

The floating center of divergence O23 is defined as the crossing point of the 
two curves C2 and C3. By definition, 

ZP1O23P2 = ZP2O23P3 = ZP3O23P4. 

Thus, the three angles spanned by the four consecutive points define a diver- 
gence angle, 

d2 = ZP2O23P3. (8) 
Distance from O23 is denoted as 

r2 = |023P2|- 

In general, the divergence angle dn and the radius r„ are evaluated from four 
consecutive points Pn-i, Pn, Pn+i and Pn+2; 

dn = LPnOnn+lPn+li (9) 
Tn = \Onn+lPl\- (lO) 
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These quantities are compared with 



= LPnOoPn+U (11) 

and 

= \Ojtl (12) 

which are defined in terms of an arbitrarily chosen center Oq. 

Given coordinates of the four points Pn {n = 1, 2, 3, 4), the position of the 
fioating center O23 is parametrically represented by means of two parameters 
X and Y- 

W3 = XiV? + YP^,. (13) 
The two parameters are determined by the two equations ([6]) and ([7]), namely 

IO23-F1P + 1023^2^ — l-Pl-f^l 1 023^2 P + IO23P3P — I-P2-P3I 



and 



i^23^T[^23^] IO23-F2I IO23-F3I 



|023-^P + 1023-^3^ — 1-^2-^31 _ |023^P + IO23P4P — I-P3-P4I 
IO23-F2I IO23-F3I 1023-^3! IO23-F4I 



(14) 



(15) 



by the cosine formula in trigonometry. 

For example, consider vertex points P„ (n = 1,2,3,4,5) of a regular 
pentacle, whose Cartesian coordinates are given by (x„, ?/„) = (cosnrf, sinnrf), 
where d = 27ia in radians and a = |. Let the center O23 of the four points Pn 
(n = 1, 2, 3, 4) be represented as (fT3l) . For this special case, it is not difficult 
to find the solution X = Y = ~ 0.2764 by geometrical considerations. 
The golden ratio 

v^+1 ...^ 
r = (16) 

is the irrational number quintessential to the phenomenon of phyllotaxis. The 
estimate of X ^ F ~ 0.3 can be used as an initial guess for the numerical 
search of solutions in general cases. 

The manner in which the fioating center fioats around may be illustrated 
by means of a line segment connecting the fioating center and the middle 
point of the middle two leaves defining the center, that is, the line segment 
O23M23, where M23 is the middle point of P2 and P3. Let us call a graph of 
the line segments a divergence diagram. If the center is fixed in space, all the 
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Figure 2: Four points Pi {i = 1,2,3,4) from a deformed Lucas pattern for which two 
curves C2 and C3 as defined in Fig. [1] fail to cross. Thus, it may happen that the center 
of divergence is not defined when divergence angle is small and variable. 



line segments radiate from the fixed center. If it is not, the hnes cross with 
each other. See Fig. Hj for instance. 

When a pattern with small divergence angles is deformed significantly, 
it may happen that no center is defined because the two curves C2 and C3 
do not cross (Fig. |2]). Even in such center can be defined formally 

by selecting four points properly. Here it is remarked only that a pattern 
consisting of more than three points may not have a definite center. In other 
words, it is very special for numerous leaves comprising a phyllotactic pattern 
to have a unique, fixed center. Assuming the fixed center is not at all a trivial 
matter. 



2.2. Floating center of parastichy 
2.2.1. Spiral phyllotaxis 

For a spiral pattern with one leave per each node, leaves are orderly 
indexed by an integer n, called the plastochron. Four points for a center need 
not be consecutive in plastochron. A divergence angle is evaluated from four 
vertex points P^, Pn+q, Pn+q' and Pn+q+q' of a unit cell of a parastichy lattice, 
where q and g' are opposite parastichy numbers. Parastichies of logarithmic 
spirals are discussed in Appendix A , where it is shown that divergence angle 
d = 2'KaQ of an ideal pattern satisfies the inequalities flA.12|) in terms of two 
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Figure 3: From four nearby points Pn, Pn+q, Pn+q' and Pn+q+q' making a tetragonal cell 
of a lattice of q and g' parastichies, a center of parasticliy oi''''' is determined by ([T7)l 
and dlH]). By definition, four angles subtended by four sides of the tetragon give a unique 

divergence angle against On''^ \ which is called the floating center of parastichy. 



auxiliary integers p and p' satisfying pq' —p'q = 1, flA.13|) . With this in mind, 

a center of parastichy 0^^'"^ •* is defined such that 

g'ZP„Oi^''^')p„H_, + qlP,,Oi^^'^"^Pn+,' = 27r, (17) 

and 

LP^o^:!^'''^p^+, = ip^^,,oir'^Pn+g+,'. (18) 

where LPnO''^''^'^ Pn+q signifies the angle subtended by line segments 0^''^'^Pn 
and O!^''^ ^Pn+q- See Fig. [3l Then, four angles made by the four points give 
a divergence angle; 

4"''?') = {27rp-lPnOi^'^'^Pn+q)/q 

= {2np- IPn+q'Oi^''^'^ Pn+q+q') /q 
= {lPnOl^'^'^Pn+q'+2np')/q' 

'iPn+qO';^''^'^ Pn+q+q' + 2 V) (19) 



The rationale behind this definition may be understood by expressing the first 
equation as IPnO^''' ^Pn+q = 27ip — dl^''^ ^q, which is to be compared with the 



7 



denominator of (lA.6p . This is the net angle between the two successive points 
along a q parastichy P„ and Pn+q when divergence angles between the two 
points are equal to d^^''^ \ The difference in the signs in front of geometrical 
angles like IPnO^^''^ ^Pn+q in f|T9|) is due to the convention that the angles 
including d^^''^ ^ are regarded as positive quantities. 

Coordinates of the center O^'"^ ^ are represented in terms of two parame- 
ters r and 6 as 

Kolr^ = r cos ^iwt + r sin OK^n- (20) 

Accordingly, f lT7|) and f lTSj) are regarded as the defining equations for r and 6. 
As is clear from Fig. [3|, it is generally assumed that r > and < ^ < |. By 
regarding 6' as a given constant, the first equation f|T7|) determines r for the 
given value of 6, or r is determined as a function of 6. Then, the parameter 
9 is fixed by the second equation ( 1T8|) . Thus, the center O!^''^ ^ for given Pn, 
Pn+q, Pn+q' 8Jid Pn+q+q' is fixcd numerically. The equations depend not only 
on the coordinates of the four points, but also on the parastichy numbers q 
and q' , while either p or p' is needed to evaluate the divergence angle in (fT9|) . 

A caveat: In the above, it is assumed that P„, with the smallest index n 
is positioned nearest to the center (Fig. [3]). In the opposite convention, P„ is 
farther from the center than Pn+q+q' ■ To adapt the above results to this case, 
the four points P„, Pn+q, Pn+q' and Pn+q+q' should be replaced by Pn+q+q', 
Pn+q', Pn+q and P„, respectively. A set of equations f|T7|) . f|T8l) . and (|T9|) is 
invariant by this replacement, whereas ( l20l) should be read as 

Pn +q+q'On''^ — T SlU 9 Pn+qPn+q+q' + "^^ COS 9 Pn+q' Pn+q+q' ■ (21) 

In the biological literature, leaves on a stem are numbered in the order of 
appearance, while primordia on an apex are counted in the opposite order. 

While dn in iQ based on four distant points is insensitive to individual 
displacements of the points, d^^'"^ •* in flTI?]) based on nearby points is insensitive 
to a collective displacement of the points. The latter has a merit of wider 
applicability in practice. It is stressed again that parastichies generally do 
not have a definite center; the center 0^^'"^ ^ thus determined depends on the 
plastochron n and the parastichy pair (g, q'). 

2.2.2. Multijugate phyllotaxis 

A multijugate pattern bears more than one leaves at each node. There 
has been no established way of systematize multijugate leaves, particularly 
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because a special index to distinguish leaves at a node is lacking. Therefore, 
in the first place, a systematic index system for a multijugate pattern is 
proposed below. Then, the center of parastichy is defined similarly as in the 
last subsection. 

A J-jugate pattern has J leaves at a node. There are J fundamental 
spirals correspondingly. Each of J leaves at the n-th node is specified by 
polar coordinates 6'„j), where the jugacy index j = 0, 1, 2, ■ ■ ■ , J — 1 is 
introduced. Any leaf may be chosen as the origin (n, j) = (0, 0) of the index 
system. Along with the plastochron n, the jugacy index j is set in order in the 
direction of the fundamental spirals. The leaf with the index (n, j) is referred 
to by the symbol P^, or simply denoted as . For instance, see Fig. [T51 for 
a trijugate system with J = 3. For convenience' sake, the value range of 
the jugacy index j is extended to all integers; the coordinates {rnj,Onj) are 
regarded as periodic in j with period of J, namely {rnj+j, Onj+j) = {rnj, Onj) 
for all n and j. Accordingly, divergence angle 

di = ^n+l,j — dn,j (22) 

is periodic in j. As implied by the notation, depends on the leaf index 
{n,j), whereas it becomes constant, d^^ = d, for an ideal pattern. 

To put it concretely, let us consider an ideal multijugate system spiraling 
in the positive direction, i.e., d > 0. A pattern with a negative angle d < 
is the mirror image of the positive counterpart. In an ideal J-jugate pattern, 
the radial coordinate is independent of the jugacy index j, 

whereas the angular coordinate is given by 

0n,, = 27C (nao + , (23) 

where the divergence angle d = l-na^^ is represented with a dimensionless 
parameter a^. Without loss of generality, < Jao < \. For multijugate 
patterns J > 1, the other quantity of interest is displacement angle between 
neighboring leaves at a node, 

A^„ = ^„,,+i-^^„,,-^. (24) 

For the ideal pattern with fl23|l . IS?^ = 0. In general, A{ fluctuates evenly 
about zero. 
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In this index system, a parastichy is specified by a set of numbers {q,j), 
increments of {n,j). Conspicuous parastichies following nearby leaves are 
shown to have j = —p, where p is an integer near Jaoq (cf. below flA.6|l ). 
Thus, the q parastichy running through a leaf -P* connects the points P^, 
Pn+% -Pn+2g) -^n+sg! ' ' iustaucc, the trijugate pattern in Fig. [13] has 

J = 3 and {q,p) = (3,1), (5,2), etc. Note that a 3-parastichy 1° — 4'^ — 

— 10^'^ — ■ ■ ■ is equivalently represented as l" — 4^ — 7^ — 10" — ■ ■ ■, as the 
jugacy superscript j is understood modulo J = 3 by the periodic extension 
prescribed above. 

Having thus prepared, the method of the parastichy center in the last 
subsection is applied. The divergence angle rf^*-*''^ ^ of a multijugate system 
is defined by the four points P^;^, P^',, Pl+g+q'; (HZD, (IIHl) and (dl 
should read 

^VpiOfe'^Op^^P + qlP^Ot^'^P^-^: = 2vr/ J, (25) 
lP^Oi^,'^')p^-P = iP^-pi^'^'^P^^,^,,, (26) 

and 




3. Effect of uniform deformation 

Consider a regular spiral pattern in which position of the n-th leaf is given 
by Cartesian coordinates {xn,yn)- The center Oq of the polar coordinate 
system (r^°\6'^°^) is set at the origin. The azimuthal angle 6'^°^ is measured 
from the x-axis. As typical cases for a phyllotactic pattern to lose track of 
the center, two types of global deformation are considered. The deformation 
is characterized by a single parameter representing strain C. 
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Uniform deformation of the first order transforms y„) as follows; 



X 



^ 1 (1 + Ci)Xn. Xn>0 
\ (1 - Ci)Xn. Xn<0 

Vn > Vn- (28) 

Uniform deformation of the second order makes 

Xn y (1 — C2)Xn, 

Vn > Vn- (29) 

The magnitude of deformation is represented by constants Ci and C2, re- 
spectively, which are assumed significantly smaller than 1. In both cases, 
the radial component r^°^ = ^x'^ + y'^ and divergence angle 6*^1^^ — O^li suffer 
modulation periodic in the azimuthal angle The periods of modula- 

tion for the first and the second order deformation are 360 and 180 degrees, 
respectively. 

The amphtude of modulation in the radius r^°^ ± Ar^°) is given by 

^ = ^1 (30) 
for the first order deformation, and 



ri°) 2 



(31) 



for the second order deformation. 

Amplitude of modulation in divergence angle, AS^\ may be estimated 
as follows. Consider an isosceles triangle with the vertex angle d^^\ the base 
of length 2L, and the height H. Then, 

L 

By the first order deformation, (0'' becomes d^'^^ ± Ac?*^") when the height 
H parallel to the x axis is modified to (1 ^ Ci)i/. By the second order 
deformation, d is modified to ± Ad^'^) as {H, L) become ((1 — C2)H, L) 
and {H, (1 — C2)L). In either case, 

d(o)±Ad(o) 

tan c^_(l^C), (32) 
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where C is Ci or C2. Accordingly, 



(33) 



In particular, 

Arf(°) ~ 38.7C de errees 
for = 137.5 degrees, and 

Ad(°) ~ 56. de grees 



(34) 



(35) 



for ci(o) = 99.5 de grees. 

The first order deformation r nay be caused by a tendency tow ard or 
against the direction of the sun ( Kumazawa and Kumazawa jSnl)). The 
second order deformation may be caused by various reasons, real or appar- 
ent. The most possible would be due to uniform compression. It applies 
also to the case where a pattern is observed from a direction oblique to the 
perpendicular direction. If the angle of inclination is denoted as if in radian 
measure, then 1 — C2 = cosip, or C2 — V'^/S. 



4. Results 
4.I. Example 

Consider a regular pattern with d = Itx jr^ (137.5 degrees) and a = 1.073 
for ([T]) and ([2]). This is an orthogonal (5,8) system corresponding to i = 5 
in ( 1A.25I) . By deforming the pattern according to (1281) with C\ = 0.05 and 
with C2 



0.05, patterns shown in Figs. 4(a) and 4(b) are obtained, 
respectively. In Fig. |U line segments radiating from the center represent the 
divergence diagram explained at the end of Sec. 12.11 The fixed center Oq for 
divergence angle and radius r^^^ is set at the center of the original pattern. 
Divergence angles d and radii r^^\ r for the pattern of Figs. 4(a) and 4(b) 



are shown in Figs. |5] and |6l respectively. All quantities show characteristic 



variations. A type of variations as exhibited by 0?^°^ in Fig. 5(a) is frequently 
met in real systems. Fig. 5(c) shows that the variations are correlated with 
azimuth 6n- In general, deformation has a larger effect on the divergence 
angles than on the radii. Variations in the divergence angle d by means of 
the floating center in Sec. 12.11 are considerably suppressed compared with 
those of A period of 180 degrees in Fig. 6(c) is the characteristic of uni- 



form compression in one direction, namely deformation of the second order. 
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• 19 

19 




Figure 4: An ideal pattern of a logarithmic spiral with the plastochron ratio a — 1.073 
is compressed (a) by (|28l) with Ci — 0.05, and (b) by (|29|) with C2 = 0.05. Leaves are 
represented with points with integer index n. Radiating lines represent 0„„+iM„„_|_i, 
where 0„„+i is the floating center of divergence for the four consecutive leaves from the 
{n — l)-th to (n + 2)-th leaves and M„„+i is the middle point of the n-th and (n + l)-th 
leaves. It is called a divergence diagram in Sec. 12.11 
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(a) 0?^°-' and d against n 



(b) and r against n 
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(c) c?^°-' and d against 



(d) Ar'^"^/r'^'^' and Ar/r against 



Figure 5: Results for the deformed pattern in Fig. 4(a) d'-'^^ and r*^"^ are divergence 



angle and radial coordinate with respect to the center of the original (undeformed) spiral 
pattern, while d and r are divergence angle and radial coordinate with respect to the 
floating center of divergence. For the horizontal axis, n is the leaf index, and 9 is the polar 
angle, or azimuth. Ar^°^ and Ar are changes in r'-'^-' and r by the uniform deformation for 
Ci = 0.05. 



14 



139 
138 
137 

136 




2 4 6 



10 12 14 16 18 20 



(a) and d against 



141 
140 
139 
138 
137 
136 
135 



60 120 180 240 300 360 




(b) r^°^ and r against 




50 100 150 200 250 300 350 



(c) d*^"' and d against 9 (d) Ar^^'/r'") and Ar/r against 9 

Figure 6: Results for the deformed pattern in Fig. |4(b)] (C2 = 0.05); cf. Fig. [SJ 
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Figure 7: Divergence angles of the original (not deformed) pattern of Fig. |4] are measured 
against a misplaced center Oq set at the middle point of leaves 1 and leaf 5. The misplace- 
ment brings about apparent wild variations in dn '' (open circles, left), which are correlated 
with azimuth 9„ (right). Divergence angle dn evaluated with the floating center of diver- 
gence (filled circles) does not depend on the wrong center Oq; the original divergence angle 
is correctly retrieved without knowing the true center. 



Divergence angles facing the compressed direction increase apparently. For 
C2 = 0.05, gives Ad ~ 1.9 degrees, which corresponds to the amplitudes 
of d^^^ in Fig. 6(c) and Ar^'^^r^^'^ in Fig. 6(d) The amplitudes of d and Ar/r 
are suppressed and enhanced compared with d^^^ and Ar^^yr^'^\ respectively. 
The factors of suppression and enhancement are about the same value ~ 1.8. 

Even if a pattern has a definite center, misidentification of the center 
causes apparently systematic variations in d^'^\ Fig. [7] is a result for the 
undistorted pattern Ci = C2 = 0, whereas the middle point of 1 and 5 is 
chosen as the nominal center Oq. The divergence d based on the floating 
center does not depend on the choice of Oq, whereas the variations in 
exhibit a period of 360 degrees when plotted against the azimuth 6. Thus, 
the misplacement of the center can be misinterpreted as the first order de- 
formation. 

These results indicate that systematic features of variations may not be 
revealed unless they are plotted against the azimuth. 



4-2. Suaeda vera 

As a typical example of real systems, ph yllotaxis of young lea ves of Suaeda 
vera is shown in Fig. 8(a) after Fig. la of iRutishauserl ( 119981 ). The actual 
size of the pattern is about 1mm in diameter. Leaf position is optically 
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(b) Divergence angles d„ and dn ^ against 
the leaf index n 
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(c) dn and dn ^ against 6*; 
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Figure 8: Results for a normal phyll otaxis of young lea.v es of Suaeda vera, whose pattern 
shown in (a) is taken after Fig. la of Rutishauseil ( 1998f) . Radiating lines in (a) represent 
a divergence diagram (cf. the caption of Fig. |4] and Sec. 12. ip . 
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read as marked in the original figure, which is a drawn figure based on a 
microtome section. A leaf is not a point. It has shape and size. In this 
case, and presumably in most past cases reported so far, position of a leaf is 
represented by the position of the main vascular bundle when it is visually 
discernible. This paper does not take up cases in which position of leaves may 
not be specified unambiguously. A fixed center Oq for in Figs. 8(b) and 
8(c) is set according to the original figure. Statistical results for divergence 
angle are = 136.8 ± 5.8° by the fixed center Oq, d = 136.9 ± 1.8° by the 
fioating center of divergence, and S^'^^ = 137.0 ± 0.9° by the fioating center 
of parastichy. It is remarkable that no significant deviation from the ideal 



limit angle of 137.5° is found, although the divergence diagram in Fig. 8(a) 



indicates that the pattern does not have a well-defined center. In the figures, 
d^' ^^ for n = 15 is ona itted because no solution was found. 

Rutishauserl (jl998[ ) has remarked angular width of leaf dbTC , 1 , clS another 
quantity of interest. From the original figure, it is estimated as i = 66 ± 13°. 
Oddly enough, the standard deviation of the angular width of leaf arc is 
larger than that of divergence angle. If the right tip of leaf arc is regarded 
as a representative point of the leaf, d^^^ = 137.4 ± 6.0°. If the left tip is 
used instead, d^'^^ = 136.1 ±5.8°. Thus, divergence angle does not depend on 
which part of leaf is regarded as a representative point. The results indicate 
that absolute positions of the leaves are not affected even though the angular 
width of each leaf arc fluctuates widely. 

Apparent lack of stability in nominal divergence angle d^^^ (open circles 
in Fig. 8(b) ) is due to systematic variations of the type expected for the flrst 
order deformation (Fig. 8(c)). The azimuth plot of d in Fig. 8(c) (closed 
circles) indicates a slight indication of peaks at ~ 80° and ~ 80 + 180°, 
signifying deformation of the second order. 

A semi- log plot of r„ and r!^^ in Fig. 8(d) conflrms the exponential rule 
Straight lines labeled with (3,5), (5,8) and (8, 13) represent the exponential 
growth according to (1A.25I) for i = 4,5 and 6, respectively. In the original 
flgure, leaves are in contact with each other along three contact parastichies 
3, 5 and 8. The growth rate log a has to be estimated in order to evaluate the 
phyllotaxis index (P.I.) in (lA.30p . This may pose a methodological problem. 
The radial growth r„ of a real system is neither continuous nor monotonic 
in n. There are systematic variations in r^°) and r„, particularly because the 
pattern does not have a deflnite center. Accordingly, a bad method or a bad 
ch oice of leaves may give an unwanted result log a < for some values of 
n. 



Rutishauserl (119981 ) has evaluated a ~ 1.06 based on two ratios for three 
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Figure 9: Results for a control shoot o f Xan thium, after the bottom photo-micrograph in 
Fig. 6 of iMaksymowych and EricksonI (jlQTTl ). Compare with Fig. [TOl 



leaves, namely 2, 18 and 31 in Fig. 8(a) (30, 14 and 1 in the original figure). 
In the presence of systematic variations, a curve fitting method objectively 



gives a more accurate result. As shown in Fig. 8(d) , the least-squares fitting 
method with the function roexp(— nloga) gives log a = 0.0580 ± 0.0014 (a = 
1.0597 ± 0.0015) and log a = 0.0593 ± 0.0023 (a = 1.0611 ± 0.0024) for r^") 
and r„, respectively. Therefore, P. I. =4. 18 by (lA.30p . The pattern is close to 
an orthogonal (5, 8) system with P.I. =4. 

4-3. Effect of gibberellic acid on Xanthium 

Figs. M and [TU] ar e results obtained after the bottona and top photomi- 
crographs in Fig. 6 of IMaksymowych and EricksonI (Il977l ). which are control 
and gibberellic acid (GA) treated vegetative shoots of Xanthium, respectively. 
The mean and standard deviation of divergence angle are rfW = 137.73±8.03° 
and d = 137.74 ± 1.46° for the former, while = 138.71 ± 19.93° and 
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Figure 10: Results for a gibberellic acid treated shoot of X anthium, after the top photo- 
micrograph in Fig. 6 of iMaksvmowvch and EricksonI ( 19771) . Compare with Fig. |9l 
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Figure 11: Lucas phyllotaxis on a capitulum of Helianthus tuberosus after Fig. 2 of 
Rvan et alJ (flQQL) . 



d = 138.55 ± 3.09° for the latter. The pattern of the control shoot in 



Fig. 9(a) indicates a well-defined center, while the treated shoot in Fig. 10 (a 



does not have a center. As a result, of the treated plant shows wild 
fluctuations. Thus, the GA treatment not only decreases the growth rate 
logo (the slope of Figs. 9(d) and 10(d)), but destabilize divergence an- 
gle appreciably (Fig. 10(b)). Nonetheless, the treatment does not affect 



the mean divergence angle. These results are consistent with analysis by 
Maksymowych and EricksonI ( 1977 ) based on the exponential growth ([2]). 
The exponential growth is corroborated qualitatively, if not quantitatively. 



as shown in Figs. 9(d) and 10(d), where reference lines labeled with paras- 
tichy pairs of Fibonacci numbers (Fj, -Fj+i) represent the exponential growth 
with the constant rate according to (IA.25|) . 

4.4- Helianthus tuberosus 
4.4-i- Lucas system: (11,18) 

In Figs. [TT] and [T2l results are shown for a capitulum of Helianthus tubero- 
sus after Fig. 2 of iRyan et al.l (Il99ll ). The size of the pattern is about 
13mm in di a meter . Here and below, digitized position of florets is taken after 
Ryan et al.l (119911 ). although the numberings are different from the original 
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Figure 1 2: Results for a capitulum of Helianthus tuberosus plotted in Fig. 1111 (After 
Fig. 2 of lRvan et all (|l99l[) V 
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figures. This is a Lucas system with the limit divergence angle of 



1 

3 + r- 



(36) 



or d 
by 



27rQ;o = 99.5°. A nominal center Oq is set at the center of gravity, or 

5:o;^ = 0. (37) 



The center of gravity is a good choice for the fixed center especially for a sys- 
tem comprising a large number of pattern units. Results for divergence angle 
are dS^^ = 99.72 ± 5.35°, d = 99.42 ± 1.28° and d^^^'^^^ = 99.47 ± 0.15°. Re- 
mark the accuracy attained in d^^^'^^\ Thus, the center of parastichy comes 
into its own in a high phyllotaxis pattern. In a closely packed pattern, young 
seeds in a central region are amenable to irregular displacements caused by 
inward pressure due to old seeds in an outer region. This is observed as wide 
fluctuations of d^^^ for n > 70 in Fig. 12(a) As shown below, the secondary 



displacement appears to be a primary source of variations in divergence angle 



of a packed pattern. The divergence angles in Fig. 12(b) indicate system- 
atic variations of the compression type. The amplitude Ad^^^ ~ 4° implies 
C2 — 0.07 by although the deformation is apparently indiscernible from 
Fig. [m In general, deviations of r„ from an exponential dependence give rise 



to shifts in parastichy numbers. In Figs. 12(c) and 12(d) , curves labeled with 
(11, 18) and (7, 11) represent the exponential dependence ([2]) with flA.lip for 
^ ~ 18' ^ ~ u f ~ f' ~ n' respectively. The limit divergence 
angle of is used for aQ. The n-dependence of r„ is neither exponential 
nor square root It is rather close to linear (Fig. 12(d)). Accord- 
ingly, conspicuous parastichies change continuously from (11, 18) near the 
rim to (7, 11) near the center. The apparent up shift of the parastichy pair 
caused by gradual decrease of log a is called rising phyllotaxis. The pattern 
in Fig, im shows falling phyllotaxis, the downshift of the parastichy pair from 
(11, 18) to (7, 11), and even down to (4, 7). It should be remarked that new 
cell primordia arise from the rim towards the center, not vice versa as often 
assumed wrongly in theoretical models (cf. @). Mathematically, the shift of 
parastichy corresponds to the fact that the function logr„ is convex upward 
as shown in Fig. 12(c) Despite the gradual shift in parastichy numbers, di- 
vergence angle is not affected (Fig. 12(a)). If phyllotaxis is to be judged by 



the divergence angle, there is no sign of change in phyllotaxis. 
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Figure 13: A rare trijugate 3(3, 5) phyllotax i s on a sectioned capitulum of Helianthus 
tuherosus. Adapted after Fig. 3 of iRvan et al. ( 1991 ). Seeds are indexed as according 
to the numbering system (n, j) explained in Sec. 12.2.21 



4.4.2. Trijugate system: 3 (3,5) 

Figs. [13] and [H] are results for a r are trijugate patte rn (J = 3) of He- 
lianthus tuherosus based on Fig. 3 of iRyan et al.l ( 119911 ). The size of the 
pattern is about 8mm in diameter. In terms of the nominal center Oq set by 
(EZD, one obtains 3(i^(°) = 137.4 ± 14.7°, 137.7 ± 12.9° and 135.3 ± 19.2° for 
j = 1,2 and 3, respectively, where 



is divergence angle depending on the jugacy index j. Note that j = 3 is 
equivalent to j = 0. In Fig. [131 the seed point is denoted as . For every 
integer n, the trijugate pattern has three organs of approximately 120 degrees 
apart from each other, which are distinguished by the jugacy index j. When 
averaged over j, nominal divergence ang le is 3rf(°) = 136.8 ± 9.0°, showing 



a large standard deviation. As shown by dH^^^ in Fig. 14(a) , variations in 
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Figure 14: Results for a capitulum of Helianthus tuberosus plotted in Fig. [T3] (after Fig. 3 
of iRvan et all (fl99lh V (a) d^^"^ and Sn^ = (d^^^^ + dl^°^ + dT^)/3, and (b) rl!°\ 



rn^^^ and r^^"^ are plotted against the leaf index n. They are evaluated in terms of the 



nominal center Oq set by ff37\ 



divergence angle of the multijugate system are substantially larger than that 
of a simple spiral system. 
The radial coordinate 



.J(0) 



\OoP^\ 



curves 
4 and 



depends linearly on n with less remarkable fluctuations. In Fig. 14(b) 
labeled with 3(3,5) and 3(5,8) are drawn by flA.261) for J = 3,z = 
J = 3, z = 5, respectively. The results suggest that fluctuations are due 
to displacements in the angular direction. Indeed, large fluctuations A^"^ = 
± 6.85° are observed for fl2^ . Multijugate systems are characterized by 
this large angular fluctuations. For this reason, the method of the center 
of divergence is barely applicable. By means of the parastichy center in 
SecEXH we get 3rf^(3'5) = 136.53 ± 2.55°, 137.85 ±2.88° and 137.43 ±3.78° 
for j = 1,2 and 3, respectively. In total, 3d^^'^^ = 137.28 ± 3.09°. The 
standard deviations are significantly reduced but still an order of magnitude 
larger than the spiral system of the last subsection. The absence of a center 
of rotation in a real J-jugate system is manifested by the very fact that 
expected rotational symmetry of order J is badly approximate. This may 
not be apparent at a glance of Fig. [131 but it is immediately checked if 
Fig. [13] is overlapped with itself after being rotated by 120 degrees, or, e.g., 
by comparing three pairs of (2°, 7^), (2^, 7^) and (2^, 7°). The rotated pattern 
does not coincide with the original one. 
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Figure 15: Dev eloping florets on a capitulum of Helianthus annuus after Fig. 4 of 



Ryan et alj (|199lD . 



4^.5. Helianthus annuus 

Figs. [15] and [16] ar e for developing fl orets on a capitulum of Helianthus 
annuus after Fig. 4 of iRyan et al.l (Il99ll ). The actual size is about 9 cm in 
diameter. Numerical results for divergence angle are = 137.49 ± 1.64°, 
d = 137.49 ± 0.56° and rf^^^'^^) = 137.50 ± 0.07°. In numerical calculation, 
some inner seeds fail to give the solution for the ce nter of para s tichy, owing 
to significant displacement from regular position. [Ryan et al.l ( 1l99ll ) have 
noted large fluctuations of open circles in Fig. 16(a) The present study 



reveals that the fluctuations are not random but mainly systematic variations 
due to local deformation of the pattern, as indicated by the azimuth plots in 



Figs. 16(b) and 16 (g 



The radius in Fig. 16(e) is fitted with 



r„ = roVl — bn + cn^ (38) 
to give ro = 0.528 ± 0.002 (arbitrary), h = (4.88 ± 0.07) x 10"^ and c = 
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Figure 16: Resu lts for a capitulum of Helianthus annuus in Fig. [15] (after Fig. 4 of 
Rvan et al.l (|l99ll )). 
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Figure 17: Displacement from regular positio n of florets on the capitulum of Helianthus 
annuus in Fig. [15] (cf. Fig. [HI). After Fig. 4 of iRvan et a"D (|l99l[ ). 



(5.01 ±0.27) X 10- 



Thus, the area per seed, 
d 



d(- 



-n] 



-(7rr^) = 7rro^(6 + 2c(-n)), 



depends on the seed index n. The positive value of the coefficient c signifies 
that inner seeds are smaller than outer seeds, as expected physically. The 
relative change in size is 2c/h = 0.002, namely 0.2 percent per seed. Accord- 
ingly, the innermost seed for n = 287 is estimated to be about half the size 
of the seeds around the rim. This is indeed seen from the original figure. 



Scatterings of divergence angles of inner seeds (Fig. 16(a)) suggest that the 



physical effect to reduce the size of seeds affects their position and divergence 
angles accordingly. In terms of the continuous monotonic function f l38|) . the 
n-dependent plastochron ratio is evaluated as 



log a = -— logr„ 
an 



b — 2cn 



2(1 — bn + en?) 



(39) 



with which the phyllotaxis index is calculated by ( ]A.30p . In Fig. 16(f) 
the divergence angle dn is shown against the phyllotaxis index. Note that 
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P.I. = 5,6 and 7 represent (8,13), (13,21) and (21,34) orthogonal systems, 
respectively. As the capitulum develops, the phyllotaxis index decreases and 
variance in d increases accordingly (falling phyllotaxis). However, the mean 
value of d is constantly held at the limit divergence angle. The deviation 



Ar„, = 



\/l — bn + on?' normalized by r„ is plotted against the az- 



imuthal angle in Fig. 16(h) Variations in r„ show a similar dependence of 
9 as dn in Fig. 16(g), namely the dependence by deformation of the second 



are com- 



order (Sec. [3]). Amplitudes of variations in Figs. 16(g) and 16(h) 
pared with Figs. 6(c) and 6(d) for C2 = 0.05. Displacements of seeds are 
explicitly shown in Fig. [T71 where each vector represents shift in position of 
a seed from its regular position determined by the mean divergence angle, 
virtually equal to the limit divergence angle. The figure indicates that (i) the 
pattern is compressed horizontally (stretched vertically) and (ii) inner seeds 
are displaced appreciably by inward pressure. These are consistent with the 
above observations, namely systematic modulation of divergence angle and 
decrease in size of inner seeds. The latter is also consistent with the obser- 



vation that seeds near the center are even aborted (iRyan et al.l (Il99ll )). 
is transformed to 



n] 



1 + 



n 



no 



(40) 



By definition, an integer is an index number of the innermost seed at the 
center, r^^ = 0. In practice, however, the index ric of the last seed may not 
be determined without ambiguity, for it depends on position of a nominal 
center and whether or not to count aborted seeds. Be that as it may, the 
numerical order of the index system can be reversed by formally replacing 
ric — n with n. Consequently, ( HOj) is equivalent to 



roWn 1 + 



n 



no 



(41) 



in the reversed index system. Fig. 16(e) shows that the radial co mponent r-n 



i s bett er fitted with ( HTl) than the square-root dependence of (jl]) . iRyan et al 
fll99lh .mve investigated products of algebraic, logarithmic, and exponential 
dependences of r„. The functional form of fj4T|) has not been considered 
previously, although it has wider practical applicability. A new constant 
no in f HTj) roughly signals the index n of a seed for which the square-root 
dependence (jl]) begins to fail; r„ oc ^/n for n <^ no whereas r„ oc n for 



29 



7 28 49^^ 36 ^ oq 

• • • • 




31 
• • • 




Figure 18: An artichoke capitulum after Fig. 9A of iHotton et al.l (|2006l ). 



n ^ Uq. As inner seeds n ~ are likely to be displaced significantly, the 
limit exponent e (~ 0.5) of the power dependence near the center r„ oc n'^ 
{n — 7- 0) would be of little physical significance. As a rough estimate, we get 
no ~ 800 for this sample, while no ~ 20 for Fig. 14(b) 



4-6. Artichoke 

Figs. [T8] and [T9] are results for an artichoke capitulum based on Fig. 9A of 
Hotton et al.l (12006). The size is about 4mm in diameter. Fig. [18] closely re- 
sembles Fig. [151 although the sizes differ by more than an order of magnitude. 
Position of primordia, which are round shaped and closely packed, is read 
from the original figure with sufficient accuracy. Numerical results obtained 
are = 137.35 ± 5.03°, d = 137.53 ± 0.52° and d^^^'^i) = 137.50 ± 0.08°, all 
falling onto the limit divergence angle of flA.22p . The variations in are 
substantially suppressed in d and d^^^'"^^^ based on the floating centers. The 
standard deviation of d^^^'"^^^ is particularly noteworthy. The nominal center 
Oq for d^^^ is determined so as to minimize the standard deviation of d^^\ 



30 





50 100 150 200 250 300 



60 120 180 240 300 360 



a) dn and dn^ against n (b) dn and dn^ against 9. 



(0) 









(34,55) 








(21,34) 






(13,21) 

















(34,55) 








(21 ,34) 






(13,21) 





50 100 150 200 250 300 350 



50 100 1 50 200 250 300 350 



(c) logr^^ and logri^'* against n (d) rn and against n 



(0) 




60 120 180 240 300 360 



(e) dn against 9 




(f) against n 




0.05 - • 



50 100 150 200 250 300 350 



(g) Ar„/r„ and Avn^ /r^^"^ (h) Displacement from regular 
against 0„ position 



Figure 19: 
(l2006l) 1. 



Results for the artichoke capitulum in Fig. [TH] (after Fig. 9A of iHotton at al 
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Still, shows relatively large variations. If the center of gravity is chosen 
for Oo, one obtains d^'^^ = 137.39 ±5.03° instead. Therefore, large deviations 
m are not due to misplacement of a fixed center. As a matter of fact. 
Fig. 19(b) clearly indicates that deviations in divergence angle are due to 
systematic variations. As a rule, apparently wild variations of divergence 
angle that cannot be suppressed by a judicious choice of a center should be 
suspected as a sign of no definite center by deformation. Fig. 19(e) reveals an 
unexpected result that remnant variations in dn are primarily not of random 
origin but of the compression type with modulation of Ad ~ 1°. As shown 
in Fig. 19(f), the radius r„ shows a square-root dependence of @, namely 
r„ = tq^/I — hn. The radius unlike the angle d^^\ is generally insensi- 
tive to the choice of a nominal center Oq, so that r,^°^ is not affected severely 
even if the nominal center Oo is chosen inappropriately. In Fig. 19(g), the 



deviation from the fit Ar„ = r„ — tq^JI — hn is plotted against the azimuthal 
angle 9. The systematic variations in r„ have a different 6'-dependence from 
that of dn in Fig. 19(e) This is inexplicable by uniform deformation consid- 



ered in Sec. [3l As shown in Fig. 19(h) , displacement of primordia from their 
regular position is not uniform; displacement vectors display local structures 
resembling vortexes of incompressible fluid. Here again, it is concluded that 
the local displacement of closely packed organs gives rise to the systematic 
variations of r„ and divergence angle dn- 



5. Discussion 



Discontent not only with untested abstractions in the mathematical lit- 
erature but also with vague expressions like "the largest available space" and 
"about 137.5°" in the botanical literature underlies this work. It is often 
assumed implicitly that a spiral pattern has a definite center. This is not 
necessarily true, as remarked in Sec. 12.11 As shown in Sec. 12.21 each cell 
of a parastichy lattice may have its own center. The center of parastichy 
is not fixed in space, but it may float around in the pattern as primordia 
arise one after another. The floating center can be the main source of ap- 
parent variations of divergence angle. In the literature, a logical leap is often 
made from a measured value of about 137.5 to the irrational number 360/r^ 
without subjecting data to statistical analysis. The lack of quantification is 
partly because the center against which to measure the angles is uncertain. 
Probably for this reason, little attention has been directed to the numerical 
accuracy of divergence angle. The standard deviation of divergence angle 
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conveys no less important information than the mean value. According to 
an observation, the standard deviation decreases sys tematically as th e plant 
grows day by day while the mean is kept constant (Williams ( 19741 ) ). This 
behavior is not confirmed by snapshot patterns investigated in this paper. 
The present work finds no significant drift of the mean divergence angle, 
which puts possible mechanisms of phyllotaxis under constraint. 

There are two cases where it is invalid to assume for a phyllotactic pattern 
to have a fixed center. In one case, the position of center is displaced as a 
leaf primordium is initiated. In the other case, leafy organs are displaced by 
secondary effects after they are initiated at their regular position. The effects 
of the former and the latter are taken into account by means of the floating 
centers in Sec. 12.11 and Sec. \2.'2\ respe ctively. It is note d that the former 



concept conforms with observation by iGreen and Baxterl (119871 ). who have 
made a general remark that the center of apical area moves as a new leaf is (or 
leaves are) initiated. Their observation is veri fi ed co ncretely by the methods 
and results of this paper. iGreen and Baxterl ( 119871 ) note that a trajectory 
of the moving center may be used as a diagnostic for phyllotactic patterns, 
which is to be contrasted with classification by me a ns of geometrical unit of 
ideal patterns ( Meicenheimer and Zagorska-Marek ( 19891 )). 

At a glance of an impressive phyllotactic pattern, our attention is apt to 
be directed to Fibonacci numbers, because our eyes tend to follow conspic- 
uous parastichies. However, parastichies do not serve much for quantitative 
purposes. Polar coordinates (r„, 9n) of pattern units cannot be deduced un- 
equivocally from parastichy numbers alone, although parastichy numbers are 
determined uniquely from coordinates (r„, 6'„). The polar coordinates are the 
most proper tool to specify a two-dimensional pattern quantitatively, only 
with which it can be investigated whether and how the divergence angle 
9n — On-i and the radial coordinate r„ depend on the leaf index n. If a 
pattern is deformed to lose a well-defined center, the meaning of divergence 
angle blurs apparently. Still, parastichies may remain httle affected owing 
to their insensitivity to quantit ative details. Even an oval- shaped capitulum 
preserves parastichy numbers ( ISzymanowska-Pulkal ( 1l994j )). The method of 



Appendix B for indexing primordia remains valid irrespective of whether 



the pattern has a center or not, because it is based solely on parastichies. 
Although a real system generally do not possess a well-defined center, the 
method of the fioating centers enables us to infer coordinates (r„. On) of the 
original pattern, as shown in the last section. 

The results of this work underline the universal characteristic of phyl- 
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Figure 20: Real and fake phyllotaxis to illustrate accurate control of divergence angle d. 
Densely packed florets arise spirally in numerical order from the rim toward the center. 
Inner florets may be omitted without affecting arrangement near the rim, that is at issue. 

(a) Normal phyllotaxis of a (21,34) pattern for d = 137.5°, quite common in nature. 

(b) Seemingly plausible but actually improbable phyllotaxis of a (21,29) pattern for d = 
136.8°. Plants distinguish the difference of less than 1° without fail, whereas the human 
eye would find it difficult to tell (a) from (b). The difference comes to the fore when the 
21st floret arises on either side of the 0th fioret; plants know the correct position (angle) 
in advance. Central portions after, say, the 120th floret are overlapped nearly perfectly by 
a rotation of about 90°, so that they are indistinguishable by means of lower parastichy 
numbers there. 



lotaxis: divergence angle during steady growth is stably held at 137.5°, or 
rarely 99.5° and other special values. This is not just a general tendency. 
As a pivotal empirical rule deduced from experimental facts that lays the 
foundation of mathematical analysis and thus underlies this intriguing phe- 
nomenon of phyllotaxis, the author believes that it is appropriate to attribute 
the status of a natural law to the universality of the accurate control of di- 
vergence angle (cf. Appendix A). To convince that this is not exaggeration, 
let us take two examples. The first is given by high phyllotaxis of florets on 
a capitulum. According to (lA.12p . a (21,34) pattern is obtained if and only 
if divergence angle d is constrained within a narrow range of ^ < 3^ < ||, 
or 137.14° < d < 137.65°. Imagine what happens if divergence angle had 
been out of the range of width of 0.5°. As a mathematical consequence, it 
is shown that (21, 29) or (13, 34) should obtain instead of the normal phyl- 
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lotaxis (21,34) (Fig. [20]) . Neither of these anomalous patterns is observed 
actually. Thus, not only the constancy of divergence angle but also the ac- 
curate control of it to the special value of 137.5° should be taken seriously as 
an empirical fact. The point is not just that a normal (21, 34) phyllotaxis is 
observed, but rather that anomalous patterns of (21, 29) and (13, 34) phyl- 
lotaxis do not coexist, or mixed, with the normal pattern. Thus, plants must 
know precisely the correct value of divergence angle beforehand. The second 
example is given by low phyllotaxis of young leaves on the apex. As shown 
in Sec. 14. 2[ divergence angle is independent of sizes of leaves. This is an 
important observation, because it apparently conflicts with the assumption 
of dynami cal models tha t diver gence angle is variable depending on the size 
of leaves (IMirabet et al.l (120121 )). So to speak, the whole is more than the 
sum of its parts. It seems ra ther appropr i ate to regard divergence angle as 



a given constant parameter (jSadoc et al.l (120121 )). Along with a review on 



various approaches to phyl lotaxis, a mod el conforming to these observations 
has been given previously (Okabe ( 20121 )). 

A remark is in order. In the ontogeny, divergence angle generally be- 
gins with 180° or 90° of a distichous or decussate pattern, and thereafter 
shows some transient fluctuations before attaining toward a universal con- 
stant value. Variations in the transient regime are substantially larger (typ- 
ically more than about 10°) than the small variations in the steady-growth 
regime investigated throughout this paper. Therefore, the two variations are 
distinguished quantitatively. The former may be understood as a natural re- 
sult of the limited numbers of apical meristem cells pushing and shoving with 
each other. The latter variations, however, are too small to be explained as 
such; actually, as remarked above, it appears rather that divergence angle is 
regulated steadily at a predetermined value, despite any adventitious factors. 
Generally speaking, living systems are so complex that large variations are by 
far easier to understand. For a quantitative understanding, the riddle of the 
small variations cannot be overemphasized. It should be investigated for its 
own sake, separately from the large transitory fluctuations during ontogeny. 

It is also a finding of this paper that apparent deviations from the univer- 
sal mean are not due to random errors as commonly expected, but they are 
actually systematic variations. On the one hand, they are likely to be due to 
directional correlation with respect to the sun (IKumazawa and Kumazawa 
(Il97ll )). For lack of space, it is only noted here that the azimuthal corre- 
lation, as plotted in Fig. 8(c), is a common phenomenon confirmed quite 
generally. On the other hand, systematic variations are caused also by lo- 
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cal displacement when organs are closely packed. Concerning the latter, 



Ryan et al.l ( 1l99ll ) have attached importance to persistent fluctuations of in- 



dividual divergence angles (open circles in Fig. 16(a)) as a real phenomenon. 



As indicated in Figs. 16(b) and 16(g), the fluctuations systematically cor- 
related with the azimuth, the absolute direction, should be regarded as a 
subsidiary phenomenon. Thus, one should be cautious not to take appar- 
ent fluctuations literally as evidence for or against a theoretical model of 
phyllotaxis. This is an important point because almost all experimental and 
theoretical divergence angles reported so far have not taken into account the 
possibility that the center for the angles is not fixed in space. Therefore, vir- 
tually all the reported results have to be reexamined carefully. If variations 
turn out to be due to misidentification of the center, they have no real sig- 
nificance. If they refiect real effects, they act to disturb the regularity. The 
disturbance, if any, conveys important information, which may be analyzed 
with a model from an appropriate theoretical perspective. 

As illustrated by Figs. [17] and 19(h), variations of divergence angle are 
interpreted naturally as necessary consequences of local distortion due to 
contact pressure between closely packed organs. The contact pressure may 
affect phyllotaxis not only quantitatively but even qualitatively. Primor- 
dia under high pressure from the surrounding primordia may fail to grow 
normally, as illustrated by innermost florets on a capitulum in Fig. [T71 In 
general, organs in a crowded portion are most likely to be aborted, or obliged 
not to arise normally. Ill effects due to an aborted primordium may prop- 
agate along a parastichy of vascular connection, along which nutrients are 
transported. The other possibility is that new primordia may happen to 
arise to flll in gaps opened between existing primordia. The author sus- 
pects that aborted or inserted primordia might ap pear as a crystall o graph ic 
defect in an otherwi s e reg ular parastichy lattice ( iZagorska-Marekl ( ll994j )). 
Szymanowska- Pulka has made an interesting observation that the 



parastichy number decreases quite more often than it increases. 

Individual divergence angles measured against a nominal center may ap- 
pear to vary very wildly. Nonetheless, it has been empirically known that 
the limit divergence angle of 137.5° is immediately obtained by evaluating 
the average of several successive angles covering a few full turns. This em- 
pirical fact is explained consistently by the observation of the present work 
that the apparent variations are not random but systematically correlated 
with the azimuth. If the variations are of random origin, the standard de- 
viation should depend on the number of sample leaves. The present study 
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has found no such dependence. In contrast, the systematic variations cor- 
related with the azimuth are averaged away within a few full turns of the 
azimuth. For the samples analyzed in this paper, the systematic variations 
are of the order of 1 degree. If this secondary effects of the directional cor- 
relation are compensated for, the accuracy and stability of angular control 
by nature should stand out all the more strikingly; it should turn out to be 
of the order of 0.1 degree as a general rule. This is an unexpected result. 
As remarked above, the numerical accuracy of this quantitative phenomena 
of living organs should arrest more attention than eye-catching Fibonacci 
integers do. 

As shown by the results, the angular equation ([T]) with constant diver- 
gence angle d remains valid quite generally. In contrast, the radial equation 
dlj) is not valid quantitatively, whereas ([2]) holds true at the shoot apex in 
a good approximation. Unlike the divergence angle d, however, the plas- 
tochro n ratio a may be changed artificially (IMaksymowych and Erickson 
( 119771 )) (Sec. 14. 3p . As a matter of fact, it has been reported that the plas- 
tochron ratio a changes naturally during shoot development. In particular, 
the ratio decreases appreciably during transition to flowering without chang- 
ing the phyllotactic sequence determined by the divergence angle d. These 
observations suggest that the growth in the radial direction r„ is controlled 
independently of the angle 9n (cf. ([3]) and ([5])). In other words, the regu- 
larity in Tn may not be dealt with on the same basis as the regularity in 6'„. 
Indeed, the present work finds no correlation between the angle 6n and the 
radius r„. This point has bearing on how to determine the numerical order 
of leaves, namely the leaf index system n. Two index systems based on the 
numerical orders of r„ and 6'„ need not be identical. Owing to the regularity 
in 6n, the index n is orderly set according to the angle 6n- In contrast, the 
radius r„ in practice is not a monotonic function of n. Fig. 8(d) indicates 
that Tn decreases non-monotonically, or the order in n (horizontal axis) is 
not preserved for r„ (vertical axis). Actually, this is normally the case when 
the plastochron ratio a is very close to 1 (high phyllotaxis), for irregular 
fluctuations in r„ outweigh regular changes by the plastochron ratio a. Even 
the he ight order of successive leaves on a stem may be interchanged ( iFujita 
fll942h l 

Excepting such minor irregularities, which are more or less expected for 
living organisms, observations generally support validity of mathematical 
description in terms of the polar coordinate system. This is not at all trivial. 
Just as an elliptic orbit of a planet pins down a special point in space, the 
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sun at the origin, the general regularities manifested plainly by means of the 
polar coordinate system signify the existence of a special singular point, the 
origin of the coordinate system. Unlike the solar system, there is no obvious 
sign at the origin of a phyllotactic pattern. Indeed this was a motivation of 
the present work; even without assuming the center, a point of the sort has 



been indicated definitely (see Fig. 9(a)). The origin or the center should be 



identified with the growing point biologically, though it is meant here in a 
narrower sense than in general botanical use. The author speculates that 
biochemical properties of the singular point, the growing center of a spiral 
pattern, might hold a key to understanding not only mechanisms of radial 
growth but also regulation mechanisms of divergence angle. 



Appendix A. Logarithmic spiral and parastichy concept revisited 

Logarithmi c spirals in p h yllota x is has bee n inve st igated iii t ensive ly by 



ChurchI fll904f) . Ivan ItersonI fllQOTl ). iRichardsl fll95lf ). lEricksonI fll983f ) and 



Jean! (119941 ). among others. The main purpose of this section is to derive the 
formulas used in the main text, flA.llI) . flA.Mp . flA.25p and f lA.26p for the 
plastochron ratio a when two opposite parastichies are orthogonal, as they 
have not been presented before in these forms. The following derivation has 
a merit of transparency by which a physical assumption and a mathematical 
approximation are distinguished by excluding unnecessary assumptions and 
complications. As remarked below, this distinction is essential to a clear 
understanding of empirical rules of phyllotaxis and Richards' phyllotaxis in- 
dex. Before that, the prevalent concept of parastichy must be made clear 
and definite. It is also an important aim of this section to point out that 
conventional expositions are incorrect or insufficient. 

In the polar coordinate system (r, 9), a logarithmic spiral through a point 
r, 9) = (1, 0) is given by 

(Al) 



where 6 is a constant. The logarithmic spiral is characterized by the property 
that the angle between the radial vector from the origin and the tangent 
vector is held constant at every point (r, 9) on the spiral. 



d9 

r— = tanyj. 
dr 



(A2) 



Hence it is also called the equiangular spiral. The angle (p is related to the 
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coefficient b in (lA.ip by 

b = cotip. (A.3) 

The fundamental spiral of phyllotaxis given by ([1]) and ([2]) forms a loga- 
rithmic spiral with 

6=1^ = ^. (A.4) 
d 27rao 

In the second equation, a reduced divergence angle ao is introduced by 

d = 2'KaQ. (A.5) 

Owing to the periodicity in angle, one may assume either 0<ao<lor— i< 
c^o ^ I without loss of generality. In the latter convention, the direction of 
the fundamental spiral is determined by the sign of divergence angle. In 
what follows, it is assumed Q < O-o < \, i.e., the fundamental spiral runs 
outward counterclockwise. It is straightforward to adapt the following results 
to spirals for — | < < 0. 

The concept of parastichy is widely used in the literature. As pointed out 
explicitly below, prevalent expositions are not satisfactory. According to a 
general consensus, a spiral connecting leaves with the indices differing by an 
integer q is referred to as a g-parastichy. However, there are infinitely many 
such spirals. Therefore, the term q-parastichy is not appropriate in the first 
place. 

As a g-parastichy, consider a logarithmic spiral r = e^^ running through 
the 0-th and g-th leaves, namely (r, 6) = (1, 0) and (a'^, 27rgao). The former is 
met from the outset, whereas the latter determines the coefficient b. Owing 
to the periodicity in angle, the polar coordinates of the g-th leaf (r^, 9q) = 
(a'',27rgao) is equivalently represented as {a'^,2'7iqao — 27ip), where p is an 
arbitrary integer. Substituting the latter into the spiral equation r = e^^ , 

loga« loga 

(A.6) 



27r (gao - p) 2-k ( «o - - 



As the spiral thus determined runs through all leaves with the index n divis- 
ible by the integer g, it is a g-parastichy. The g-parastichy depends on the 
integer p. There are as many g-parastichies as the number of integers. The 
fundamental spiral with flA.4p is a special case of flA.6P for (p, g) = (0, 1). In 



the literature, the most conspicuous, shortest one is usually regarded as the 
g-parastichy. Therefore, the integer p is identified with the integer nearest to 
gao- As a matter of fact, this need not be the case. 
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In practice, it is necessary and sufficient to let p be an integer approxi- 
mating a number qa^. The denominator of the middle expression of (1A.6P 
represents the net angle between two successive leaves on the g-parastichy. 
To keep the angle within a reduced range of width 27r, a multiple of a full 
turn, 27rp, is subtracted from the nominal angle of 27rgao- Hence the integer 
p is the winding number of the fundamental spiral executed between consecu- 
tive leaves on the g-parastichy. If the number qa^ is not an integer, p is most 
properly chosen to be either the largest previous or the smallest following 
integer to qao. In the former case, qa^—p > 0, the q parastichy spirals in the 
same direction as the fundamental spiral. In the latter case, qao —p<0, the 
direction of the q parastichy is opposite to the fundamental spiral. If qa^ is 
an integer, then p = qa, and ip = by { \A.3\i . In this special case, the paras- 
tichy is not a spiral curve but a straight half-line radiating from the center. 
Therefore it is especially called the orthostichy. It goes without saying that 
orthostichy in the strict sense does not exist in real life, as it does not make 
sense to ask whether qa^ — p is zero or not; it is equivalent to asking whether 
ao is an irrational or a rational number. In any case, the parastichy depends 
not only on the parastichy number q but also on the winding number p, or 
more precisely on ao ~ | whose sign and magnitude uniquely determine the 
sense and slope angle of the parastichy. 

As the parastichy, or (lA.6p . depends only on the fraction ^, it is sufficient 
to consider the irreducible pair of p and q. That is to say, if p = 0, then 
q = 1; otherwise, p and q are restricted to coprime integers, i.e., p and q 
are not evenly divided by any integer greater than 1. By the assumption 
< ao < ^, it is sufficient to consider irreducible fractions between and |. 
For the denominator g up to 8, they are 

p_011112112313 

g ~ T' 2' 3' 4' 5' 5' 6' 7' 7' 7' 8' 8' ^ 

Hence, parastichies for g = 5, 7, 8 are not unique. By way of illustration, 
five parastichies for | = | and | are shown in Fig. IA.21I for ao = It 
is not common to notice 5-parastichies in the pattern of Fig. IA.21I This is 
because the pattern has more conspicuous parastichies. A numerical measure 
of conspicuousness is given by the absolute value of ao — I, which is a winding 
angle per leaf of the parastichy. In a word, the less winding the parastichy 
is, the more conspicuous it appears to the eye. In other words, the better 
the fraction | approximates the divergence angle a^, the more conspicuous 
the parastichy is. To find conspicuous parastichies for a given value of ao. 
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Figure A. 21: A phyllotactic pattern for ao = jq and a = 1.1. Two distinct 5-parastichies 
connect the same set of points, 1, 6, 11, 16, 21; the solid curve is a 5-parastichy for 
{p,q) — (1,5), and the dashed curve is another 5-parastichy for {p,q) — (2,5). The latter 
spiral winds in the direction opposite to the fundamental spiral. 



a sequence of irreducible fractions arranged in the numerical order, called a 
Farey sequence, is of much help. By arranging (lA.Tp . we obtain the Farey 
sequence of order q = 8; 

011111213231 

T'8'7'6' 5'4' 7'3' 8'5' 7'2' ^ 

As ao = ^ hes between | and | in this sequence, it is properly understood 
that 3 and 7 parastichies make a conspicuous pair for the pattern of Fig. lA.2ll 
In spite of the above remark, the conventional term of the g-parastichy is 
used throughout this paper in order not to complicate matters by introducing 
a new term like a ^g-parastichy. In practice, we get along without encounter- 
ing the ambiguity because the number qao always happens to be nearly an 
integer, namely a Fibonacci number, by the empirical fact that divergence 
angle is given by a special irrational number (cf. ( IA.22|) ). Therefore, the 



integer is always identified with p unwittingly. Another reason is that dif- 
ferent sets of parastichies with a common parastichy number never become 
conspicuous at the same time, because adjacent terms in a Farey sequence 
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have different denominators ( Hardy and Wrightl (jl979l )). This is a conse- 
quence of (1A.13P and concretely checked by ( lA.Sp . 

The second shortest g-parastichy m ay become rele vant. For instance, 
consider for ao = and q = 3 (ISwintonI (|2012| )). The 3-parastichy 

with p = is actually not a 3-parastichy because (p, q) = (0, 3) is reducible 
to (0,1), the fundamental spiral. Therefore, a parastichy for (p, g) = (1,3) 
should rather be referred to as the 3-parastichy. Thi s example r e futes the 
preval ent statement that p is the integer nearest to g«o ( Jean ( 19941 ): SwintonI 
(l2012h l 

Along with the q parastichy with b in (]A.6|) . let us consider a g'-parastichy 
for another arbitrary integer q'. In a similar manner as above, a g'-parastichy 
is given by a logarithmic spiral with 



b' 



logo 



27r (^ao 



(A.9) 



for b in (lA.ll) . where p' is an integer approximating g'ao- According to 
(IA.3p . the coefficient b and b' for the g and g' parastichy are related with 
the slope angle ip and ip' hj b = cot ^9 and b' = cotip', respectively. When 
the g-parastichy and the g'-parastichy are mutually orthogonal, the identity 
if + if' = ^ holds true. Accordingly, bb' = cot ip cot ip' = —1, i.e.. 



log a 



log a 



271 {ao - f ) 271 (ao 



-1, 



or 



log a = 2n 



\ 




(A.IO) 



(A.ll) 



which is real and positive. The inequality bb' < 0, which holds when the 
two parastichies are opposite, signifies that ao should lie in between the two 
fractions | and ^7. Let the former fraction be numerically larger than the 
latter, without loss of generality. Then, 



P ^ 

— <ao< -. 
q' q 



(A.12) 



In the case of ao < 0, both p and p' become not positive, so that the in- 
equalities in (IA.12P should be reversed, whereas (lA.lip remains valid formally 
as it is. 
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Figure A. 22: A phyllotactic pattern with divergence angle ao = t^^ and the plastochron 
ratio a in (jA.lip for ^ = ^ and ^ — \- This pattern has an orthogonal parastichy 
pair of (<?, g') — (3,13), i.e., 3-parastichies (dotted) and 13-parastichies (dashed) cross 
orthogonally. Nonetheless, it is usually referred to as a (8, 13) phyllotaxis, because 8- 
parastichies (solid) with ^ = | are as conspicuous as 13-parastichies. As pq" — p"q = 1 
for (g", q) = (8, 13) whereas pq' — p'q — 2 for (q', q) — (3, 13), the former has a reason to 
be regarded as the fundamental pair. Parastichies are imaginary constructs, though some 
may have real observable effects. 
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A pair of parastichies is commonly referred to by the denominator pair 
(g, q'). The parastichy pair {q, q') is usually not chosen so that the parastichies 
are nearly orthogonal; it is chosen to be irreducible, i.e., such that the two 
limiting fractions in (IA.12|) constitute successive terms of a Farey s e quenc e. 
See Fig. IA.22I By a theorem of number theory ( iHardy and Wrightl ( 119791 )). 
it is equivalent to stating that integer pairs for a pair of parastichies {p, q) 
and (p', q') are chosen to satisfy 



' ' 1 
pq — p q = 1. 



(A.13) 



The identity is confirmed for adjacent fractions in f lA.8|) . The irreducible pair 
of parastichies should not be confused with the irreducible numbers {p, q) for 
a g-parastichy discussed above. The former may be rephrased as conspicuous, 
although the term carrying ordinary connotations might better be avoided. 
In Fig. IA.22t the parastichy pair (3,13) is not conspicuous (irreducible), 
whereas not only (8, 13) but (3, 5), (5, 8), (13, 21) and so on are conspicuous 
(irreducible). Thus, (IA.13|) should be considered as a convenient prescription 
for good choices of the parastichy pair. Physically, the left-hand side of f lA.13p 
represents the number of leaf points in a unit cell of the lattice spanned by 
q and q' parastichies (Fig. IA.22|) . f lA.13|) means not only that q and q' are 
coprime integers but that p and p' are determined independently of uq. As 
a result, the parastichy numbers do not depend on the divergence angle, 
insofar as the condition (IA.12P is met. Note that the phyllotactic pattern 
is completely characterized by the two parameters, the plastochron ratio a 
and the divergence angle ao, which are not determined uniquely by parastichy 
numbers alone. The relation between the divergence angle ao and the limiting 
fractions and ^ to be used for flA.lip is immediately read from a diagram 



presented in lOkabd (120 111 . 120121 ). 

An opposite parastichy pair of a multijugate system has the greatest 
common divisor J > 1, so that the pair is written J(g, g'), where q and q' are 
coprime integers. Instead of flA.lip . one obtains 



log a 



27r 



p 



— JaQ JaQ — 



p 

7 



(A. 14) 



where p and p' are integers near Jaoq and Ja^q' . The divergence angle cto 
of the J-jugate system satisfies 



P ^ J ^P 
q' q 



(A.15) 
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If J(g, q') is a conspicuous pair, 



Jq'p — Jqp = J, 



(A.16) 



or flA.131) holds true. Thus, the divergence angle ao and the plastochron ratio 
a of a multijugate system may be better expressed in the forms of Joq and 
a"' to compare them with a single spiral system (Sec. I4.4.2p . 

In the above, no approximation has been made besides the basic equations 
and ([2]). Above all, the derivation is not based on any ideal assumption 
on leaf shape, e.g., contiguous circles abundant in the phyllotaxis literature. 
This is practically of crucial importance because any argument based on 'per- 
fectly circular leaves' must meet with severe criticism from experimentalists. 
As a matter of fact, parastichy is a property of a point lattice, so that it is 
independent of shape and size of pattern units. If flA.lSP is to be re garded as 
the condition for a visible parastichy pair, this is what iJeanI (119941 ) calls the 
fundamental th eorem of phyllo taxis. But then, one must be careful about 
what is visible (Swinton ( 20121 )). See a remark below (IA.13P on the term 
irreducible. 

The above formulation is general. According to number theory, linear 
Diophantine equations, (1A.13|1 and ( IA.16|) . are solved for any parastichy pair 
of integers. If the parastichies are orthogonal, the plastochron ratio a is given 
by (lA.lip or (]A.14|) . Nevertheless, nature adopts only selected numbers as 
the parastichy pair. Parastichy numbers of a normal phyllotaxis comprise 
the Fibonacci sequence. 



for 



Fi — 1, 1, 2, 3, 5, 



^= 1,2, 3,4,5, 6,---, 



(A.17) 



respectively. The Fibonacci numbers satisfy the mathematical identity 

- i^„2F,+i = (-1)\ (A.18) 
Therefore, a solution of (1A.13P is given by 

{q,q',p,p') = (F,+i,F„i^_i,F,_2) (A.19) 
if i is an even integer, or by 

{q,q',p,p') = (F„F,+i,F,_2,F,_i) (A.20) 
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if i is odd. In either case, (lA.lip gives 



log a = 271^ 



1-2 



(A.21) 



for the orthogonal Fibonacci parastichy system of f lA.2ip is valid 

insofar as lies between and tt^- Nonetheless, it is usually taken for 
granted that divergence angle ao is identified with the mathematical limit of 



^ and 



"i+i 



- as i — )■ oo, that is. 



/ + r 

flRichardd fll95lh : Ijeani fll983l )). where r is the golden ratio given in (fT6|) . 
The divergence angle of c? ~ 137.5° for (]A.22p is called the limit divergence 
angle of the main (Fibonacci) sequence. To show that ( ]A.22l) is the limit of 
the fraction is expanded with respect to after the exact formula 



i-r)- 



V5 



is substituted for the denominator and numerator; 



_i-2 



i-r) 



-i+2 



_ (-1)^757 



-2i 



(A.23) 



(A.24) 



The second term vanishes in the limit i — t- oo. Substituting 
(1X241) into flA:2Tl) . 

logo 



22]) and 



(A.25) 



for the orthogonal Fibonacci system (Fj, -Fj+i) with the limit divergence 
angle (IA.22p . This result is generalized to an orthogonal J-jugate system 
J(Fj,Fj_|_i) with the limit divergence angle = t^'^/J, 



log a 



In terms of common logarithms, (1A.25P is expressed as 

277^5 



logio (logio «) = -^logio(^ ) + log 



10 



^rlog(lO)^ 



(A.26) 



(A.27) 
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which is transformed into 



log,o(-^) " logio(r^) ■ ^ ^ ^ 

Expressed in numbers, 

i - 1 = 0.37918 - 2.39249 logio (logio «) • (A.29) 



The right-hand side is the phyllotaxis index (P.I.) of Richards f 195ll ). Thus, 



it is proved that the phyllotaxis index is nothing but the index i of the 
Fibonacci system (Fj, -Fj+i) minus one. This is an integer by definition. In 
natural logarithms, the index is given by 

^ log(2vrv/5/r) _ log(loga) _ ^ 
log log 

P.I. may take any value if it is regarded as a function of log a. In applying 
the formula flA.30P to real systems, it should be kept in mind that log a in 



the logarithm must be a positive number, namely a > 1 by assumption. 
The above derivation clearly indicates that Richards' phyllotaxis index is 



based on the physical assumption ( ]A.22I) and the math ematical a pprox ima- 



tion ( [Al24l) . This is not obvious from the derivations by iRichardsl (119511 ) and 



JeanI (119831 ). Interestingly, it is mostly the latter approximation that turns 



out to be less reliable, i.e., the rational number for small i cannot be 



replaced by the irrational number r ^. The former assumption (lA.22p . or 



d ~ 137.5°, is empirically supported quite accurately, as shown in the main 
text. In the author's view, the regularity expressed by ( IA.22I) . or On = nd 
with d = ±27r/r^, should be called the fundamental law of phyllotaxis and 
treated particularly as such. Let us incidentally remark that logarithmic spi- 
rals of phyllotaxis are far more miraculous than similar spirals appearing in 
certain growing forms like nautilus shells by the very fact that the angle d is 
not only held constant but specifically fixed at 27r/r^ = 137.5°. 

For a non-logarithmic spiral pattern, parastichy numbers will be shifted 
in a Fibonacci sequence depending on part of the pattern. For instance, the 
square-root spiral by (jl]) is formally regarded as having the plastochron ratio 
a depending on the leaf index n, namely logo = ^(n-i) ^ogn by r„ = a"~^ = 
y/n, when a reference scale is set as ri = 1. 
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Leaves on a stem may be represented in terms of a cylinder coordinate 
system as 



0„ 



nh, 
nd. 



(A.31) 



where h is the length of an internode relative to the girth and d is a divergence 
angle. Accordingly, the fundamental helix is parametrically given by 



27r 



-9. 



where 



h_ 

"0 



(A.32) 



(A.33) 



It is easily shown that the angle (f that the helix (]A.32p makes with the z axis 
is expressed in terms of b in (]A.33I) by the same equation as (lA.3p . namely 
b = cot (f, Therefore, one may follow the same derivation as (1A.11|) to obtain 



h 



\ 




p 

7 



(A.34) 



for the orthogonal parastichy system (g, q'). Note that the plastochron ratio 
a of a spiral pattern on the apex and the internode length /i of a cylindrical 
pattern on the stem are formally related by 



2ti 



log a. 



(A.35) 



Although the cylindrical representation is frequently used in mathematical 
studies, real systems do not obey (lA.3ip because the internode /i, unlike d 
and a, is not constant even approximately. ( 1A.34P should be understood as 
a reference result of an abstract model. 



Appendix B. How to number sunflower seeds 

To measure divergence angle, the sequential order of organs has to be 
identified according to their age or plastochron. Sometimes this may in- 
volve laborious tasks, particularly for a high phyllotaxis pattern like packed 
seeds on a sunflower head. In principle, all seeds reachable from a reference 
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/ 


12 3 4 


mod(/i?, Q) 


0/^2/'4\l/3\ 


g = 5 


2 2 12 1 



Table B.l: Numbers to make a primordia front circle (see Fig. IB.231) . This table is for a 
parastichy pair of {Q,Q') = (5,8). In the middle, mod(/i?, Q) denotes the remainder of 
IR divided by Q, where R — 2Q — Q' — 2 and an integer / runs from to Q — 1 = 4. The 
last row is a 1-2 sequence for Q = 5, which is obtained by positing either 2 or 1 depending 
on whether the remainder increases or decreases. 



seed along parastichies are numbered by successively adding or subtracting 
parastichy numbers. A drawback of this simple method is that a single mis- 
calculation spoils all the following numbering. Therefore, it is practically of 
much help to have a systematic device. The following method starts from 
choosing a layer of seeds that lie near a rim circle, which are regarded as a 
"siege" to start numbering whole seeds inward. In what follows, seeds are 
counted from the outermost toward the center. 

Let us consider a pattern with opposed parastichies {Q,Q'), where in- 
tegers Q and Q' have no common divisor. For definiteness, Q and Q' are 
assumed to satisfy 1 < Q' /Q < 2, which practically holds true in most cases. 
(By this assumption, capital letters are used for the parastichy numbers.) 
The bounding integer 2 plays a key role below. In addition, an auxiliary 
integer R = 2Q — Q' is introduced. For every integer / = 0, 1, 2, ■ ■ ■ Q — 1, 
calculate the remainder of the division of IR by Q, which is denoted as 
mod{IR, Q). See Table ED for (Q, Q') = (5, 8) and R = 2Q-Q' = 2. In the 
last line, either 1 or 2 is inserted between columns depending on whether the 
remainder mod(/-R, Q) increases or decreases as / increases by 1. The last 1 
at the right bottom is set because mod(2J, 5)= for J = on the left end is 
smaller than 3 on the right end for 1 = 4. There are R ones and Q — R twos, 
so that the numbers add up to R + 2{Q — R) = Q'. Thus, the 1-2 sequence 
in the last line of Table IB. II represents a partition of Q'{= 8) into Q{= 5) 
parts (2+2+1+2+1=8). 

With this table at hand. Fig. IB.231 explains how to fix initial seeds on a 
front circle from which to start numbering the inner seeds. The sequence of 
ones and twos in the last line of Table iBTTl represents the number of steps that 
have to be made in the direction of Q parastichies before shifting back one 
step in the minus direction of Q' parastichies. Starting from a reference seed 
numbered 0, the seed 2 is reached after two steps of Q and minus one step of 
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Q' . Then it goes to the seed 4 after the same steps, because 2Q — Q' = +2. 
The next is the seed 1 by one step along Q and one step back along Q' 
(4 + Q — Q' = 1). In this way, the front circle is closed as it returns to the 
seed after five shifts of —Q'. The front circle thus defined wit hout reference 



to a n ominal center does not agree with the primordia front of iHotton et al 



( 120061 ) , which is defined such that a new primordium added to this front lies 
closer to the center than all primordia belonging to the front. 

Cartesian coordinates of seeds may be harves t ed fr o m a d i gital p hoto 



image by means of a digitizing software flMitchelll (l2009h : iKiiskI fl201ll )l In 



so doing, useful formulas are derived if one decides to pick a constant number 
of seeds from each of Q parastichies. Let this number be denoted as X^ax- 
While collecting data in a one- dimensional sequence, the starting position in 
each Q parastichy has to be successively shifted according to the 1-2 sequence, 
as described in Fig. IB. 231 Every seed in the sequential data is assigned 
two-dimensional coordinates {X,Y) set along the {Q,Q') parastichies, where 
X = 0, 1, 2, ■ ■ ■ , Xjnax — 1 and F = 0, 1, 2, ■ ■ ■ , Q — 1. In the one-dimensional 
sequence, the order of the seed at (X, Y) is given by 

N = X + X^^,Y, (B.l) 

while the plastochron of the seed is 

Plastochron Index : n = XQ + mod{RY, Q), (B.2) 

which is not to be confused with the phyllotaxis index in (lA.30p . (IB.ip is 
inverted to give 

X = mod(X,X^ax), 

Y = int(iV/X^ax), (B.3) 

where int(iV/Xmax) = {N — mod(X, Xmax))/Xmax means the largest integer 
that does not exceed X/Xmax- By substituting (]B.3|) into (IB. 20 . the index 



number n of the X-th seed is obtained without manual calculation. In (]B.2p 



the initial seed at (X, Y) = {0, 0) is set to have n = 0. The initial seed is ar- 
bitrary but it should be chosen properly so as to make the front circle as large 
as possible while keeping the circle within the capitulum. Seeds remaining 
outside the front are indexed with negative numbers. Once completed, the 
index system may be renumbered at one's discretion. 

The above method of data collection makes use of the 1-2 sequence. If one 
were interested not in numbering, but only in picking seeds to make a front 
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4 



Figure B.23: How to make an encircling layer of seeds on a parastichy lattice with a 
parastichy pair of (Q, Q') = (5, 8). See Table IRD and the text. 



51 



circle, then it is sufficient to know a periodic 1-2 sequence. The periodic 
1-2 sequence has period Q; all sequences 22121, 12122, 12212, 21221 and 
21212 for Q = 5 are regarded as equivalent. A periodic 1-2 sequence may be 
deduced recursively by simple rules to replace 2 with 21 and 1 with 2; 

1 ^ 2. 

2 ^ 21, (B.4) 
For the main sequence of phyllotaxis, starting from 

Q = 3: 2 2 1, (B.5) 

one gets 

5 : 21 21 2, 



then 

8: 21 2 21 2 21, 

and so on. The rules in (lB.4p may remind Fibonacci's original problem, 
in which 1 and 2 correspond to new and mature pair of rabbits. The 1-2 
sequence in the correct order is obtained by reversing the sequence before 
applying flB.4|) . It goes as follows: 



3 : 221 

122 
5 : 22121 
12122 
8 : 22122121 
12122122 
13 : 2212212122121 
1212212122122 
21 : 221221212212212122121 
121221212212212122122 
34 : 2212212122122121221212212212122121 
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(a) 



29 



1650 I 
3 37 71 
24 5892 
114579 

1^74 
40 95 
6 61 
27 82 
48 

1436890 
1 



84221 34 13 
63'?76 5|g684726 5 
[ 97 81 60 39 



94 



73 



52 
86 31 
65 10 
99 44 
78 23 



90 77 98 

43 64 85 ^,93 nu^fi- 
9 30 51^72^^59254612 



75' 

101 885420 
,80.F33 



9662 
;41 7 



,2 
,15 



(b) 



Figure B.24: Digitizing an image from I JeanI (|1994| ) by way of illustration. See the text 
for a procedure of making the encircling front layer. The direction of the fundamental 



spiral is opposite to that of Fig. 20(a) 



The 1-2 sequence depends on Q' too. For instance, it is 21211 for {Q,Q') = 
(5,7), while 22121 for {Q,Q') = (5,8). Given a 1-2 sequence for a pair, 
1-2 sequences for other pairs in the same Fibonacci sequence are routinely 
derived, according to the rule of reversal and replacement. From 21211 for 
(Q, Q') = (5, 7), the rule gives 2221221 for the next pair (Q, Q') = (7, 12) of 
the (5, 7) sequence, 5, 7, 12, 19, • • ■. The result is confirmed by the method 
of Table EH 

Returning to the subject, let us take the cover image of IJearu ( Il994j ) as 
an example. The first thing to do is to detect a shell layer of a well-defined 



encircling lattice spanned by opposed parastichies. In Fig. 24(a) , there are 
Q = 34 main parastichies involuting clockwise and Q' = 55 steeper opposite 
parastichies {R = 13). The ordered 1-2 sequence for Q = 34 is given just 
above. For simplicity, let us choose X^ax = 3 for the thickness of the front 
layer. A procedure is: (i) import the image into the digitizer, (ii) Set a 
Cartesian coordinate system with proper scales of x and y axes, (iii) Fix 
the outermost layer of the primordia front according to the 1-2 sequence. 
Check if the front closes properly (the top and bottom of Fig. IB.23p . (iv) 
Start clicking the seeds in sequence (the middle of Fig. IB.23|) . (v) Index the 
digitized Cartesian coordinates by (IB. 21) and (IB. 31) . 

The parastichy pair easiest to follow with the eye is to be used as {Q, Q'), 
but the choice is not unique. The same index system is obtained if a lower 
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parastichy pair is used, namely {Q' — Q,Q), {2Q — Q', Q' — Q), etc. 

Parastichy numbers of a multijugate system have a common divisor J, 
so that they are expressed as J{Q, Q') in terms of coprime integers Q and 
Q'. To close the primordia front of the J-jugate system, Q steps of a 1-2 
sequence are repeated J times. Meanwhile, the Y coordinate runs through 
JQ integers from to JQ — 1. The plastochron index is given by flB.2p . The 
jugacy index of the numbering system proposed in Sec. 12.2.21 shifts orderly 
along parastichies. For the seed at {X, Y), it is given by 

Jugacy Index : j = mod(-PX + P'Y, J) (B.6) 

by means of integers P and P' satisfying PQ' — P'Q = ±1. 

References 

Church, A. H., 1904. On the Relation of Phyllotaxis to Mechanical Laws. 
On the Relation of Phyllotaxis to Mechanical Laws. Williams & Norgate, 
London. 

Erickson, R. O., 1983. The geometry of phyllotaxis. In: Dale, J., Milthorpe, 
F. (Eds.), The Growth and functioning of leaves: proceedings of a sympo- 
sium held prior to the thirteenth International Botanical Congress at the 
University of Sydney, 18-20 August 1981. Cambridge University Press, pp. 
53-88. 

Fujita, T., 1942. Zur Kenntnis der Organstellungen im Pflanzenreich. Jpn. 
J. Bot. 12, 1-55. 

Green, P. B., Baxter, D. R., 1987. Phyllotacticpatterns: Characterization by 
geometrical activity at the formative region. Journal of Theoretical Biology 
128, 387-395. 

Hardy, G., Wright, E., 1979. An Introduction to the Theory of Numbers. 
Oxford Science Publications. Clarendon Press. 

Hotton, S., 2003. Finding the center of a phyllotactic pattern. Journal of 
Theoretical Biology 225 (1), 15 - 32. 

Hotton, S., Johnson, V., Wilbarger, J., Zwieniecki, K., Atela, P., Cole, C, 
Dumais, J., 2006. The possible and the actual in phyllotaxis: Bridging the 
gap between empirical observations and iterative models. Journal of Plant 
Growth Regulation 25, 313-323. 



54 



Jean, R. V., 1983. AUometric relations in plant growth. Journal of Mathe- 
matical Biology 18, 189-200. 

Jean, R. V., 1994. Phyllotaxis: A Systemic Study in Plant Morphogenesis. 
Cambridge Univ. Press, Cambridge, New York. 

Kiisk, v., 2011. Digitizer, http://www.physic.ut.ee/ kiisk/files-eng.htm. 

Kumazawa, M., Kumazawa, M., 1971. Periodic variations of the divergence 
angle, internode length and leaf shape, revealed by correlogram analysis. 
Phytomorphology 21, 376-389. 

Maksymowych, R., Erickson, R. O., 1977. Phyllotactic change induced by 
gibberellic acid in xanthium shoot apices. American Journal of Botany 64, 
33-44. 

Matkowski, A., Karwowski, R., Zagorska-Marek, B., 1998. Two algorithms of 
determining the middle point of the shoot apex by surrounding organ pri- 
mordia positions and their usage for computer measurements of divergence 
angles. Acta Soc. Bot. Pol. 67, 151-159. 

Meicenheimer, R. D., 1986. Role of parenchyma in Linum usitatissimum leaf 
trace patterns. American Journal of Botany 73, 1649-1664. 

Meicenheimer, R. D., Zagorska-Marek, B., 1989. Consideration of the ge- 
ometry of the phyllotaxic triangular unit and discontinuous phyllotactic 
transitions. Journal of Theoretical Biology 139, 359-368. 

Mirabet, V., Besnard, F., Vernoux, T., Boudaoud, A., 2012. Noise and ro- 
bustness in phyllotaxis. PLoS Comput Biol 8, el002389. 

Mitchell, M., 2009. Engauge Digitizer - Digitizing software, 
http: / / digitizer.sourceforge.net/. 

Okabe, T., 2011. Physical phenomenology of phyllotaxis. Journal of Theo- 
retical Biology 280, 63-75. 

Okabe, T., 2012. Vascular phyllotaxis transition and an evolutionary mech- 
anism of phyllotaxis. http://arxiv.org/abs/1207.2838. 

Richards, F. J., 1951. Phyllotaxis: Its quantitative expression and relation 
to growth in the apex. Philos. Trans. R. Soc. B 225, 509-564. 



55 



Ridley, J. N., 1982. Packing efficiency in sunflower heads. Math. Biosci. 58, 
129-139. 

Rivier, N., OcceUi, R., Pantaloni, J., Lissowski, A., 1984. Structure of Benard 
convection cells, phyllotaxis and crystallography in cylindrical symmetry. 
J. Phys. (Paris) 45, 49-63. 

Rutishauser, R., 1998. Plastochrone ratio and leaf arc as parameters of 
a quantitative phyllotaxis analysis in vascular plants. In: Jean, R. V., 
Barabe, D. (Eds.), Symmetry in plants. World Scientific, pp. 171-212. 

Ryan, G. W., Rouse, J., Bursill, L., 1991. Quantitative analysis of sunflower 
seed packing. Journal of Theoretical Biology 147 (3), 303-328. 

Sadoc, J.-F., Rivier, N., Charvohn, J., 2012. Phyllotaxis: anon-conventional 
crystalline solution to packing efficiency in situations with radial symme- 
try. Acta Crystallographica Section A 68 (4), 470-483. 

Swinton, J., 2012. The fundamental theorem of phyllotaxis revisited. 
http://arxiv.org/abs/1201.1641v2. 

Szymanowska-Pulka, J., 1994. Phyllotactic patterns in capitula of Carlina 
acauhs L. Acta Soc. Bot. Pol., 229-245. 

van Iterson, G., 1907. Mathematische und Mikroskopisch-Anatomische Stu- 
dien iiber Blattstellungen. G. Fischer, Jena. 

Vogel, H., 1979. A better way to construct the sunflower head. Mathematical 
Biosciences 44 (3-4), 179-189. 

Williams, R., 1974. The shoot apex and leaf growth: a study in quantitative 
biology. Cambridge University Press. 

Zagorska-Marek, B., 1994. Phyllotaxic diversity in Magnolia flowers. Acta 
Societatis Botanicorum Poloniae 63, 117-137. 



56 



